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We calculate the low-lying eigenvalues and eigenvectors of the hermitian domain wall Dirac operator on various 
gauge backgrounds by Ritz minimization. The mass dependence of these eigenvalues is studied to extract the 
physical 4 dimensional A, whose spectral density is related to ('i/'i/') through the Banks-Casher relation, and 5m, 
which represents the effects of the residual chiral symmetry breaking in domain wall formalism on a per eigenmode 
basis. The topological structure of the underlying gauge field is examined by measuring the Fs matrix elements 
between the low-lying eigenmodes. 



1. DIAGONALIZATION METHOD 

We use the conjugate-gradient method pro- 
posed by Kaukreuter and Simma |jl[|^ to cal- 
culate the lowest 19 eigenvalues and eigenvec- 
tors of the hermitian domain wall Dirac operator 
Dh = where D is the domain wall Dirac 

operator and R is the reflection operator in s di- 
rection (look at 1^ for a detailed description of 
the conventions). 

Table 1 

The bare quark masses used for the Wilson 
ensembles are 0.0, 0.0025, 0.005, 0.0075, 0.01. 
For the Iwasaki ensemble, only the lowest mass, 
0.0005, is different. 
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For each of the configurations listed in Table |l]. , 
we calculate the eigenvalue spectrum at five dif- 
ferent masses. There are two main reasons for 
this practice. First, by studying the mass de- 
pendence as described in section 2, we can re- 
move from the eigenvalue spectrum a portion of 
the residual chiral symmetry breaking effects. For 
the Iwasaki action at (3 — 2.6, the residual mass 
is very small, which makes the convergence of the 
Ritz minimization prohibitively slow at zero bare 
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quark mass. Thus, we had to use nonzero bare 
quark masses (5 x 10^^ and up) to speed up the 
convergence, and then extrapolate to the region 
of interest. 

2. SPECTRAL DENSITY AND CHIRAL 
SYMMETRY BREAKING 

In the continuum limit, the spectral density of 
Dirac operator is related to the chiral condensate 
by the Banks-Casher relation 



12V m 12V \ ^ A2 + m2 



and 



lim lim = T7;P(^)- 
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(2) 



Since Dh is a continuous function of m/ , we can 
expand and parameterize its eigenvalue Aha as 

= ril/Xj + {mj + 5m,f) + 0{m)). (3) 

From 
dDu 



druf 



(4) 



where —rajRQ^'^^ is the mass term in the fermion 
matrix, we can derive 
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The quantity ("0^) is defined on the boundary for 
domain wall fermions, and given by 




By comparing Eq. |^ and Eq. ^ we can recogonize 
the parameter as an eigenvalue of the four di- 
mensional Dirac operator. Here, Srrii represents 
the contribution to the eigenvalue from the chi- 
ral symmetry breaking effects of coupling of the 
domain walls. 




Figure 2. A and Sm distribution for Wilson en- 
semble at (3 = 6.0, Ls = 8. 
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Figure 1. A and Sm distribution for the Wilson 
ensemble at /3 = 6.0, Lg — 16. The horizontal 
line in the upper panel is p(0) calculated from 
(i/'tA) using Eq. ^. The vertical line in the up- 
per panel is the minimum of the largest A in each 
configuration, which means there is no error in- 
troduced by dropping higher modes to the left of 
this line. The vertial line in the lower panel is the 
residual mass determined independently from the 
midpoint Ward/Takahashi identity term. 

Fig. |l|, H and || are the distribution of A and 6m 
for the ensembles listed in Table |l|. A rough agree- 
ment between the peak of the Sm distribution and 
mres from the mid-point term is observed. A nice 
agreement between the spectral density and (tpi^) 
is confirmed except for the peak at A = 0, which 
is explained by the abundance of zero modes at 
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Figure 3. A and Sm distribution for Iwasaki en- 
semble at (3 — 2.6, Ls = 16. 

finite volume, and is expected to disappear in the 
infinite volume limit. The comparison between 
Fig |l] and ^ shows the expected Lg dependence 
of Sm. Fig and have the same Lg, volume 
and strength of coupling but different types of ac- 
tion, which makes the Sm an order of magnitude 
smaller for the latter. 

3. GAUGE FIELD TOPOLOGICAL 
STRUCTURE 

In the continuum limit, the Dirac operator has 
75 symmetry, which ensures that the zero modes 
have chirality of ±1, and the non-zero modes, 
with chirality 0, are paired by the application of 
75 . The number of zero modes should correspond 



3 



to the winding number of the gauge field. 






Figure 4. Lego plots of three sequential configu- 
rations separated by 3 and 4 heat bath sweeps. 

In the domain wall formalism, we display the 
matrix elements |(A//_i|r5|A/fj)| as a 3 dimen- 
sional lego plot. This provides a way of visualiz- 
ing the chiral properties of the domain wall Dirac 
operator and the topological structure of the un- 
derlying gauge field. Fig ^ shows a set of such 
lego plots for the Wilson ensemble at /3 = 6.0, 



La = 16. In the first lego plot, the lowest 4 of 
the eigenvectors are also eigenvectors, while 
the rest of them are paired. Fig || also shows that 
during the heat bath evolution, when a configura- 
tion transforms from the topological sector with 
topology 4 to that with topology 2, some complex 
configurations result in between which can not be 
unambiguously categorized as belonging to any 
regular topological sector. Tabic ^ lists the per- 
centage of this kind of complex configuration for 
each ensemble in Table ^ The readers are ref- 
ered to C. Dawson's poster for a more quantative 
study of the complex configurations. 

Table 2 

List of the values of rrires and percentage of com- 
plex configurations. 
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4. CONCLUSIONS 

Our prescription of extracting A and 5m distri- 
bution from mass dependence of the eigenvalues 
of Dh is consistent with the Banks-Casher rela- 
tion and the usual ways of determing (V'V') ^-nd 
rrires- The good chiral properties of domain wall 
fermions makes them a very good tool for study- 
ing the topological structure of the background 
gauge field. We find that this structure is closely 
related to the type of action used in generating 
these configurations. 
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